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CSFIT:  A  FORTRAN  Program 
for  Charge-Sheet  Model  Fitting  of  MOSFET  Data 

L.  C.  Witte 

Semiconductor  Electronics  Division 
National  Bureau  of  Standards 
Gaithersburg,  MD  20890 

ABSTRACT 

A  FORTRAN  program,  CSFIT,  has  been  developed  for  fitting  an  expression 
for  the  current-voltage  (I-V)  characteristics  of  a  long-channel  MOSFET  to 
experimental  I-V  curves.  The  one-dimensional  charge-sheet  model  developed 
by  Brews  provides  the  basis  for  the  I-V  characteristics.  The  I-V  characteris¬ 
tics  given  by  this  model  are  optimized  with  respect  to  a  set  of  experimental 
data  using  the  fiatband  voltage  and  the  mobility  as  the  only  adjustable  pa¬ 
rameters.  The  program  is  written  so  that  multiple  sets  of  I-V  data  can  be 
fit  simultaneously  if  desired.  The  user  must  supply,  in  specified  formats,  a 
current-voltage  data  file,  a  device  parameter  file,  and  a  starting  value  file. 

KEY  WORDS:  charge-sheet  model;  fiatband  voltage;  I-V  characteristics; 
mobility;  MOSFET;  parameter  extraction 

1.  INTRODUCTION 

A  FORTRAN  program,  CSFIT,  has  been  developed  for  fitting  an  expression  for  the 
cxirrent-voltage  (I-V)  characteristics  of  a  long-channel  metal-oxide-semiconductor  field- 
effect  transistor  (MOSFET]  to  experimental  I-V  curves.  The  sign  convention  adopted  is 
appropriate  for  n-channel  enhancement-mode  devices.  The  one-dimensional  charge-sheet 
model  developed  by  Brews  [1|  provides  the  basis  for  the  I-V  characteristics.  The  I-V  char¬ 
acteristics  given  by  this  model  are  optimized  with  respect  to  a  set  of  experimental  data 
using  the  fiatband  voltage  and  the  mobility  as  the  only  adjustable  parameters. 

The  program  CSFIT  executes  the  optimization  according  to  the  least  squares  criterion  us¬ 
ing  a  subroutine  STEP  IT,  developed  by  J.F.  Chandler  [2].  STEPIT  finds  a  local  minimuTn 
of  any  given  smooth  function  of  several  parameters.  CSFIT  reads  in  dc  MOSFET  I-V  data 
and  physical  device  parameters  and  then  calls  the  subroutine  STSET  for  the  initialization 
of  the  STEPIT  input  parameters.  STEPIT  performs  the  minimization  of  the  specified 
current  model.  The  subroutine  MODEL,  called  by  STEPIT,  uses  the  current  equation 
derived  by  Brews  [l].  The  MODEL  routine  then  calls  the  subroutine  CSO.  The  CSO  and 
FUNCTION  subroutines,  developed  by  C.  L.  Wilson  [3],  solve  for  the  potential  along  the 
oxide-silicon  surface  at  the  drain  and  soiirce  end  of  the  channel  needed  to  determine  the 
drain  current.  The  subroutine  FUNCTION  is  called  by  CSO  to  determine  the  actual  drain 
current  once  the  potential  at  the  drain  and  source  has  been  determined. 


1 


In  section  2,  a  brief  discussion  of  the  formulation  of  the  charge-sheet  model  is  presented.  In 
particular,  the  current  equation  is  given  along  with  a  disctission  of  the  method  of  solving 
for  the  boundary  conditions  of  the  potential  along  the  oxide-silicon  surface  for  a  given  gate 
voltage  and  drain  voltage.  The  current  equation  and  the  boimdary  conditions  are  precisely 
what  the  CSO  computer  code  solves. 

In  section  3,  the  specific  procedures  needed  to  implement  CSFIT  are  discussed.  The 
function  of  each  of  the  principal  subroutines  is  outlined.  Limitations  of  the  program  are 
discussed. 

In  section  4,  verification  of  CSFIT  is  presented.  An  example  of  usage  is  demonstrated 
by  optimizing  the  fiatband  voltage  and  mobility  when  radiation-induced  interface  efiects 
are  present.  The  optimized  values  are  then  used  to  generate  current-voltage  curves  in  the 
subthreshold,  linear,  and  saturated  regions  of  MOSFET  operation.  The  input  files  used 
to  obtain  the  optimized  parameters  and  the  output  file  results  are  given  in  Appendix  B. 
A  print-out  of  the  CSFIT  computer  program  is  given  in  Appendix  C. 

2.  THEORY 

The  charge-sheet  model  of  the  MOSFET  assumes  that  the  carrier  density  in  the  inversion 
layer  is  contained  in  a  charge  sheet  of  zero  thickness.  Thus,  the  current  is  constrained  to 
fiow  along  the  oxide-silicon  interface  and  there  is  no  voltage  drop  across  the  inversion  layer. 
For  long-channel  devices^  it  can  be  assumed  that  Poisson’s  equation  must  be  solved  in  one 
dimension  only  becatise  the  potential  variations  along  the  channel  are  gradual  compared  to 
those  normal  to  the  channel.  The  discontinuity  condition  across  the  oxide-silicon  interface 
is  then  applied,  implying  that  the  carrier  density  per  unit  area  N(y)  is  contained  in  a  charge 
sheet  of  zero  thickness.  The  coordinate  system  used  measures  'V  along  the  channel  of  a 
MOSFET  from  the  source  toward  the  drain,  “z”  along  the  width  of  the  drawn  gate,  and 
'be”  perpendicular  to  the  plane  of  the  gate. 

Using  an  arbitrary  potential  along  the  oxide-silicon  surface  ^(x  =  0,y),  Poisson’s  equation 
for  the  potential  in  the  oxide  and  the  semiconductor  can  be  solved.  Then,  with  a  derived 
equation  for  N(y)  (refer  to  Brews  for  details  [l]),  an  integral  equation  for  the  surface 
potential  is  obtained.  This  equation  is  then  solved,  with  the  resulting  drain  current  given 
by 

/  =  (I)  {Co*  (1  +  0Vg)  {4>sl  -  M  -  ^Co*  -  <t>lo) 

-  qNA^B  -  1)^^*  -  {P4>S0  - 

+  qN^LayM  [{P<t>SL  -  1)'/"  -  [0<t>so  -  1)'/’] }, 


2 


where  <ffsL  =  ^s{L)  and  <^so  =  ^5(0).  This  is  eq  (11)  in  the  paper  by  Brews  [l].  Definitions 
of  the  parameters  used  in  eq  (11)  are  given  in  Appendix  A.  The  value  of  ^50  and  <^sl  for 
a  given  gate  voltage  Va  and  dr^  voltage  Vd  must  be  known  before  the  current  eq  (11) 
can  be  solved. 

At  the  source  end  of  the  channel,  <^so  >>  tak^  from  the  one-dimensional  Poisson  equation 
for  zero  drain  voltage  with  (0<i>so  —  1)  >  exp{—3<^so),  yielding 

.  \3  -1 1/3 

CosiVa  -  <f>so)  =  qNALBy/2  3<f>so  “  ^  ^  j  ’  (12) 


which  is  eq  (12)  in  Brews’  paper.  Equation  (12)  can  be  rearranged  to  provide  an  algorithm 
for  <I>SQ.  The  iterative  form  becomes 

4o‘  =  +  ifn{  [^(Va  -  4o)]  ’  /  [l  +  Wso  -  1)/  =q>  ^(4o  -  Ws)] },  (13) 

where  is  the  fth  trial  value  for  ^50  and  the  initial  trial  value  is 

4is„  =  14a  -  gfa 


Once  the  iterative  algorithm  for  the  potential  at  the  source  has  converged,  the  resulting 
potential  can  be  used  in  the  current  equation. 

To  determine  the  potential  at  the  drain  end  of  the  channel,  for  a  given  Vd  and  Vq,  eqs 
(14),  (16),  and  (17)  from  the  paper  by  Brews  [l]  were  used.  These  equations  are  as  follows: 


/  \*  11/3 

qN{o)  =  qNALBy/2{  3<f>so  -  1  +  j  -  {04>so  -  1)'/'}  (14) 


0^SL  =  ^^50  +  + 


qN{L)  =  CoxiVa  -  <^sl)  —  qNAl>B'^{3*^SL  -  1)^'^^. 


Using  eqs  (14)  and  (17)  in  (16)  results  in  a  value  for  the  potential  at  the  drain  end  of  the 
channel.  With  the  source  and  drain  potentials,  ^50  and  ^3 it  the  current  eq  (11)  can  be 
solved.  Equation  (11)  is  the  function  used  to  be  minimized  in  the  optimization  routine. 
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3.  COMPUTER  IMPLEMENTATION 
A.  SOFTWARE  MODULES 


In  order  to  run  tlie  program,  the  user  should  carry  out  the  following  steps.  The  user  needs 
to  obtain  a  copyrighted  copy  of  STEPIT  from  the  distributing  organization*  and  create 
three  input  files  in  the  required  format.  After  creating  the  three  input  files,  the  user  will 
need  to  link  CSFIT  to  STEPIT.  Once  the  routines  are  linked,  the  user  may  then  proceed 
to  execute  the  optimization  routine.  The  program  will  then  prompt  the  user  for  the  I-V 
data  file  name,  the  device  parameter  file  name,  the  initial  starting  value  file  name,  and 
finally,  the  desired  output  result  file  name.  This  sample  set  of  linkage  control  statements 
has  been  used  to  execute  the  CSFIT  program  on  the  DEC  VAX  11/780  under  the  VMS 
system. 

The  computer  program,  CSFIT,  consists  of  three  sections:  a  data  input  and  initialization 
section,  a  fitting  program,  STEPIT,  developed  by  J.P.  Chandler  |2],  and  a  model  section, 
as  shown  in  figure  1.  In  the  following  paragraphs  of  this  section,  a  circled  number  refers 
to  an  encircled  step  on  the  fiow-diagram  given  in  figure  1. 

The  input  section  of  CSFIT  consists  of  the  input  routine  0  and  STSET  0.  The  input 
routine  reads  in  measured  dc  transistor  data  amd  known  device  parameters.  This  routine 
normalizes  device  parameters.  The  STSET  subroutine,  developed  by  Chandler  [2],  reads 
in  constraints  on  the  minimization  process  and  sets  controlUng  parameters  for  the  fitting 
routine  STEPIT. 

The  fitting  section  uses  STEPIT  0  to  minimis  the  objective  function  subject  to  the 
constr^ts  imposed  by  the  user.  The  method  used  in  STEPIT  is  an  accelerated  version 
of  the  ‘cyclic  variation’  or  ‘one-variable-at-a-time’  method.  First,  each  fitting  parameter 
x(j)  is  varied  individually.  If  progress  is  made  by  finding  a  lower  value  of  the  objective 
function,  the  step  size  is  doubled.  However,  the  niimber  of  steps  in  this  phase  is  limited. 
When  a  local  minimum  has  been  foimd,  quadratic  interpolation  is  used  in  an  attempt  to 
refine  the  position  of  the  minimum.  The  minimum  point  is  never  shifted  unless  a  better 
(lower)  value  of  the  objective  function  is  found. 

Once  each  of  the  x(j)  has  been  varied,  attempts  are  made  to  move  along  the  resultant  of  all 
of  these  individual  displacements.  Doubling  of  the  step  size  follows  each  successful  step, 
this  time  with  no  limit  on  the  number  of  steps.  Quadratic  interpolation  again  follows. 
STEPIT  consists  of  several  subprograms  which  are  referenced  in  [2]. 

The  model  section  consists  of  the  subroutines  MODEL  0,  CSO  0,  and  FUNCTION  0. 
The  CSO  subroutine  solves  the  charge-sheet  equations  for  the  potentials  at  the  drain  and 
source  end  of  the  channel,  passing  the  values  to  FUNCTION  which  calculates  the  drain-to- 

STEPIT  is  available  from  the  Quantum  Chemistry  Program  Exchange,  Indiana  University, 
Bloomington,  IN  47401. 


4 


.V.V-V, 


cso 

@1 

FUNCTION 

Model 

Equations 


Figure  1.  CSFIT  computer  program  sections  and  subroutine  flow  diagram. 


source  current.  Tlie  calculated  current  is  passed  via  0  and  0  to  the  MODEL  subroutine, 
which  generates  an  objective  function  which  is  the  sum  of  the  squares  of  the  differences 
between  the  measured  and  charge-sheet  model  drain  currents.  The  STEPIT  subroutine 
receives  the  objective  function  via  ®  and  proceeds  to  ®  to  determine  if  the  convergence 
tolerance  on  minimization  is  satisfied.  When  the  convergence  tolerance  is  satisfied,  CSFIT 
is  finished,  ff  the  convergence  tolerance  is  not  satisfied,  steps  0  to  0  are  repeated. 


The  first  file  is  the  current-voltage  (IV)  data  file.  The  data  are  read  into  an  (npts)  x  3 
array,  where  npts  is  the  total  number  of  points  in  all  of  the  data  sets.  CSFIT,  without 
any  modifications,  assumes  that  the  measured  source-to-drain  voltage  for  each  data  set  is 
in  npts  X  1,  the  measured  drain  current  is  in  npts  x  2,  and  the  measured  gate  voltage  is  in 
npts  X  3.  This  specific  format  can  easily  be  changed  to  meet  the  user’s  needs  if  so  desired. 
Thus,  the  required  I-V  data  file  would  be  as  follows: 


npts  =  the  total  number  of  data  points 


vds(i) 

vg(i) 

id(i) 

vds(j) 

vg(j) 

idO) 

vds(k) 

vg(k) 

id{k) 

The  second  file  is  the  device  parameter  file.  The  following  format  is  required: 


nvds 

=  the  number  of  data  sets 

vds(l) 

=  source-to-drain  voltages  for  each  data  set 

vds(i) 

i=l,nvds 

np8et(l) 

=  number  of  data  points  in  each  set 

npset(i) 

i=l,nvds 

nch 

=  average  channel  doping  (cm^) 

tox 

=  oxide  thickness  (cm) 

Ijun 

=  lateral  source  and  drain  diffusion  under  the  gate  (cm) 

Ich 

=  drawn  length  of  polysilicon  (cm) 

z 

=  width  of  channel  (cm) 

The  third  file  is  the  initial  starting  value  file,  with  the  following  required  format: 

initial  fbv(l) 
initial  mob(l) 
initial  fbv(2) 
initial  mob(2) 
initial  fbv(nvds) 
initial  mob(nvds) 

The  initial  values  of  fiatband  voltage  (fbv(i))  must  be  equsd  to  or  less  than  zero.  The  sub¬ 
routine  STSET  is  an  initialization  subroutine  for  the  actual  optimization  routine  STEPIT. 
STSET  allows  the  user  to  set  the  following  controlling  parameters: 

x(i)  =  initial  values  of  the  fitting  parameter,  on  which  the  fxmction  to  be 

minimized  depends. 
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xaiax(i)  =  the  upper  limit  on  the  value  of  x(i). 

xaiin(i)  =  the  lower  limit  on  the  value  of  x(i). 

deltx(i)  =  the  initial  step  by  which  x(i)  is  varied. 

delmn(i)  =  the  lower  limit  (convergence  tolerance)  on  the  step  size  for  x(i). 

mask(i)  =  set  nonzero  if  x(i)  is  not  to  be  varied,  but  rather  is  to  be  held  at  its 

initial  value. 

ntrac  =  determines  the  amount  of  printout  from  STEPIT. 

ntrac=l  to  obtain  trace  during  the  minimization  process 
ntrac=0  for  initial  and  file  output  only 
ntrac=-l  for  no  output  except  error  messages 

STSET  will  set  all  the  above  to  default  value  except  ntrac;  i.e.,  ntrac  must  be  set  by 
the  user.  Without  modification,  CSFIT  is  written  with  the  above  parameters  set  to  the 
following  values: 

x(i)  =  initial  starting  value  file 

xmax(i)  =  default  value  of  +1.0E+35 

xmin(i)  =  default  value  of  -1.0E>35 

deltx(i)  =  default  value  of  x(i) /lO 

delmn(i)  =  l.OE^ 

mask(i)  =  default  value  of  0 

ntrac  =  1 

The  user  need  not  be  concerned  with  modification  of  the  controlling  parameters  unless 
a  fitting  parameter  is  desired  to  be  held  constant  (mask(i)=l),  or  a  fitting  parameter  is 
desired  to  be  held  within  a  smsdler  range  of  values. 

B.  LIMITATIONS 

Often,  the  error  subroutines  in  STEPIT  cannot  be  computed.  The  error  routines  sometimes 
start  with  poor  initial  starting  values  for  the  CSO  algorithm,  resulting  in  the  square  root 
of  a  negative  number.  If  this  occurs,  the  computed  results  should  be  considered  valid. 

If  the  CSFIT  routine  is  given  poor  initial  starting  values  for  the  fitting  parameters,  the 
CSO  algorithm  will  result  in  a  square  root  of  a  negative  number,  and  hence,  the  user  must 
correct  initial  starting  values. 

4.  EXAMPLE  EXECUTION 

In  this  section,  CSFIT  is  demonstrated  using  results  obtained  with  an  irradiated  n-channel 
enhancement-mode  MOSFET.  Prior  to  irradiation  and  after  irradiation  exposure  to  10,  50, 
100,  and  500  krad(Si),  the  current- voltage  characteristics  of  the  tramsistor  were  measured. 
All  devices  were  biased  with  +9  V  on  the  gate  and  source,  drain,  and  substrate  grotmded. 
The  parameters  used  in  the  simulation  are  given  in  Table  1.  Techniques  for  determining 
these  parameters  have  been  disciissed  in  the  literature  [4]. 


Table  1 

Device  Parameters 


Lch 

6.42/im 

Vg,Vd 

0.5  to  5.0  V,  0.25  V 

tos 

0.056/im 

rich 

1  X  10i®/cm3 

Ljun 

0.70/xm 

z 

16.96^m 

Fitted  values  for  the  flatband  voltage  and  mobility  were  obtained  simultaneotisly  for  the 
0,  10,  50,  100,  and  500  krad(Si)  doses,  and  then  were  used  in  the  charge-sheet  model  to 
generate  model  I-V  curves. 

Comparison  of  the  charge-sheet  model  fitted  values  with  measured  data  shows  good  agree¬ 
ment  in  the  subthreshold  region,  the  linear  region,  and  the  saturation  region.  Fig^e  2 
compares  the  fit  to  meastired  data  for  an  unirradiated  sample.  In  figure  3,  the  fit  after  a 
10  krad(Si)  dose  is  compared  to  measured  data.  In  figure  4,  the  fit  after  a  50  krad(Si]  dose 
is  compared  to  measured  data.  As  the  dose  is  increased  above  50  krad(Si),  the  effects  of 
the  interface  trapped  charge  energy  distribution  become  important.  As  shown  in  figures  5 
and  6,  for  devices  irradiated  at  100  krad(Si)  and  500  krad(Si),  respectively,  the  accuracy 
of  the  model  degrades  with  increased  irradiation  dose. 

The  linear  region  I-V  characteristics  generated  &om  the  charge-sheet  model  and  the  fitted 
values  of  mobility  and  fiatband  voltage  are  compared  to  measured  data  as  a  function  of 
dose  in  figure  7.  Linear  region  results  obttuned  are  in  good  agreement  with  meEisured 
results. 

Charge-sheet  model  I-V  cxirves  in  the  saturation  region  using  the  fitted  parameters  are 
compared  to  measured  data  in  figxire  8  for  an  unirradiated  sample.  Figure  9  gives  a  similar 
plot  after  a  100  krad(Si)  dose. 

5.  CONCLUSIONS 

A  fitting  program  has  been  developed  to  provide  an  optimum  set  of  parameter  values 
for  fiatband  voltage  and  mobility  for  a  long-channel  MOSFET  based  on  the  charge-sheet 
model  formulated  by  Brews  [l].  This  fitting  procedure  has  been  demonstrated  using  data 
&om  an  n-channel  enhancement-mode  MOSFET  exposed  to  ionizing  radiation. 
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Figrire  2.  CSFIT  results  compared  to  measured  data  for  an  unirradiated  transistor 
Vi=0.25  V. 
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Figure  5.  CSFIT  results  compared  to  measured  data  for  a  100  krad(Si)  dose,  Vrf=0.25  V 
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Figure  S.  Charge-sheet  model  characteristics  in  the  saturated  region  using  the  fitted  pa¬ 
rameter  compared  to  measured  data  for  an  unirradiated  transistor. 


Figure  9.  Charge-sheet  model  characteristics  in  the  saturated  region  using  the  fitted  pa¬ 
rameter  compared  to  measured  data  for  a  100  krad(Si)  dose. 
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Appendix  A:  Definition  of  Parameters 
^SL  =  4s{L)  =  the  potential  at  the  drain  end  of  the  channel 
^50  =  (0)  =  the  potential  at  the  source  end  of  the  channel 

nch  =  average  channel  doping 

m  =  intrinsic  carrier  concentration 
4b  =  pin  =  bulk  Fermi  potential 

Lb  =  (Debye  Length)  = 

Co*  —  oxide  capacitance  per  unit  area 

tox  =  oxide  thickness 

Lch  ~  channel  length  ~  Lpoiy 

Lpoiy  =  gate  poly  length 

Ljun  =  lateral  Junction  depth 

Z  =  channel  width 

fjt*  —  effective  mobility  (fitting  parameter) 

VpB  =  fiatband  voltage  (fitting  parameter) 

Va  =  gate  voltage 


Vd  =  drain  voltage 


APPENDIX  B:  Example  Execution  Files 


The  following  three  required  files,  (1)  I-V  data  file,  (2)  parameter  file,  and  (3)  initial  starting 
vsJue  file,  were  used  to  obtain  the  verification  results,  with  the  (4)  restilting  output  file 
included.  The  I-V  data  file  used  consisted  of  403  points.  There  was  a  total  of  five  data 
seta,  one  set  for  each  radiation  dose.  The  files  were  as  follows: 


(1)  I-V  data  file 
403 


0.27147 

8.597500E-01 

1.299160E-08 

0.25198 

7.079650E-01 

1.428790E-09 

0.25165 

6.080440E-01 

• 

1.564180E-00 

0.25048 

■ 

4.591630E-01 

1.422954E-00 

0.25079 

9.345300£M)1 

3.364700E-0d 

Vj  V-  Id 
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(2)  parameter  file 

5  i 

0.10  Vrf(t) 

0.25 

0.5 

1.0 

2.0 

64  mpts(i) 

67 

69 

82 

121 

l.OE+16  556E-08  0.7E-04  6.416E-04  16.956E-04  rich  Ux  Ljun  Lpdy  Z 

(3)  initial  starting  value  file 

0.1  500  Vft{i)  n{i) 

0.1  500 

0.1  500 


(4)  The  resulting  output  file 

Vfb  -  -7.59070E-01  Mu  « 
Vg  Id 


0.7000E+00 
0.7500E+00 
0.8000E-I-00 
0.8500E-I-00 
0.9000E+00 
O.OSOOE-t-OO 
O.lOOOE-t-01 
0 .  lOSOE-t-Ol 
O.llOOE+01 
0 . 1150E+01 
0.1200E+01 
0.1250E-I-01 
0.1300E+01 
0 . 1350E+01 
0 . 1400E+01 
0 . 14S0E+01 
0 . lSOOE+01 
0.1550E+01 
0 . 1600E+01 
0.1660E+01 
0 . 1700E+01 
0.1750E+01 
0,1800E+01 
0.18S0E+01 
0.1900E+01 
0.19S0E+01 
0.2000E+01 
0.20S0E+01 
0.2100E-I-01 
0.2150E-»-01 
0.2200E+01 
0.2250E-t-01 
0.2300E+01 
0.2350E-t-01 
0.2400E-t-01 
0 . 2450E+01 
0.2500E+01 
0.2550E+01 
0.2600E-t-01 
0 . 2660E+01 
0.2700E+01 
0.2750E+01 
0.2800E+01 
0.2850E-)-01 
0 . 29008+01 
0 . 2950E+01 
0.3000E+01 
0.3050E+01 
0.3100E+01 
0.3150E+01 
0.3200E+01 
0.3250E+01 
0.3300E+01 
0.3350E+01 
0.3400E+01 
0.3450E+01 
0.3600E+01 


0.2205E-08 
0.5492E-08 
0.1168E-07 
0.3415E-07 
0.7210E-07 
0 . 1960E>06 
0 . 4269E-06 
0.8419E-06 
0 . 1389E-05 
0.2331E-05 
0.3363E-05 
0.4518E-05 
0.5754E-05 
0.7038E-0S 
0 . 8341E-05 
0.9659E-05 
0 . 1098E-04 
0 . 1229E-04 
0 . 1359E-04 
0 . 1489E-04 
0 . 1617E-04 
0 . 1744E-04 
0 . 1857B-04 
0 . 1994E-04 
0.2116E-04 
0.2237B-04 
0 . 2357E-04 
0.2476E-04 
0 . 2593E-04 
0 . 2709E-04 
0 . 2823E-04 
0 . 2936E-04 
0.3048E-04 
0.3159E-04 
0.3268E-04 
0.3377E-04 
0.3484E-04 
0 . 3590E-04 
0.3695E-04 
0.37998-04 
0. 39028-04 
0 . 40048-04 
0 . 41058-04 
0 . 42058-04 
0 . 43048-04 
0.44038-04 
0 . 45008-04 
0 . 45978-04 
0 . 46938-04 
0 . 47878-04 
0 . 48828-04 
0 . 49758-04 
0 . 50688-04 
0.51608-04 
0 . 52518-04 
0 . 53428-04 
0 . 54328-04 


448.10  Vd  »  0.05  Set  1 
Ideal 

0.82658-09 
0.30888-08 
0 . 10688-07 
0.30648-07 
0.87578-07 
0 . 23458-06 
0 . 55948-06 
0.11638-05 
0 . 14158-05 
0.20918-05 
0 . 28458-05 
0.36558-05 
0.45048-05 
0 . 53828-05 
0.62818-05 
0.71988-05 
0.81288-05 
0.90708-05 
0 . 10028-04 
0 . 10988-04 
0.11958-04 
0 . 12928-04 
0 . 13898-04 
0 . 14878-04 
0. 15868-04 
0.16858-04 
0.17848-04 
0.18838-04 
0.19838-04 
0 . 20838-04 
0.21838-04 
0.22838-04 
0 . 23848-04 
0 . 24858-04 
0 . 25868-04 
0 . 26878-04 
0.27888-04 
0 . 28898-04 
0.29918-04 
0.30928-04 
0.31948-04 
0.32968-04 
0.33988-04 
0.35008-04 
0 . 36028-04 
0.37048-04 
0.38068-04 
0.39088-04 
0.40118-04 
0.41138-04 
0.42168-04 
0 . 43188-04 
0.44218-04 
0 . 45238-04 
0.46268-04 
0 . 47298-04 
0 . 48328-04 


o.aesoB-i'Oi 

0.3600B+01 

0.3660B+01 

0.3700B+01 

0.3750B+01 

0.3800B-t-01 

0.3850B+01 

0.3900B-»-01 

0.3060B+01 

0.4000B-t>01 

0.4050B+01 

0.4100B-I-01 

0.4150B-»-01 

0.4200B-»-01 

0.4250B-t-01 

0.4300B+01 

0.4360B+01 

0.4400B+01 

0.4450B-t-01 

0.4S00B+01 

0.4S50B+01 

0.4600B-t-01 

0.4650E-)>01 

0 . 4700B+01 

0 . 4750B+01 

0 . 4800B+01 

0.4850E+01 

0.4900B+01 

0.4050B-I-01 

O.SOOOB+01 

0.5050B+01 

0.5100B+01 

0.6150B+01 

0 . 5200E+01 

0 . 5250B-«-01 

0.S300E+01 

0 . 5350B<f01 

0.5400E+01 

0 . 5450E+01 

0 . 5600E+01 

0 . 5650B+01 

0 . 5600E+01 

0 . S650E+01 

0 . 5700B+01 

0.5750E-H01 

0.5800E+01 

0 . 5850E+01 

0.5000E+01 

0.5950E+01 

0.6000E+01 

0 . 6050E+01 

0 . 8100E+01 

0.6150E+01 

0 . 6200E+01 

0 . 6250E+01 

0 . 6300B+01 

0 . 6350E+01 

0.6400E+01 

0.6450B-t-01 

0 . 8500E+01 


0 . 5521E-04 
0 . 5609E-04 
0 . 5897E-04 
0 . 5785B-04 
0 . 5872B-04 
0 . 5957B-04 
0 . 8043E-04 
0 . 8128E-04 
0.8212E-04 
0 . 8296E-04 
0 . 8379E-04 
0 . 8461E-04 
0 . 8543E-04 
0 . 8625E-04 
0 . 8706E-04 
0 . 8788E-04 
0 . 8867B-04 
0 . 8948E-04 
0 . 702SE-04 
0.7104E-04 
0.7182B-04 
0 . 7259E-04 
0 . 7337E-04 
0.7414E-04 
0 . 7490E-04 
0 . 7565E-04 
0.7841B-04 
0.7718E-04 
0.7790E-04 
0.7864E-04 
0.7938B-04 
0 . 8012E-04 
0.8085E-04 
0.8157E-04 
0 . 8229E-04 
0.8301E-04 
0 . 8372E-04 
0 . 8443E-04 
0.8514E-04 
0 . 8584E-04 
0 . 8654E-04 
0 . 8723E-04 
0 . 8792E-04 
0.8861E-04 
0 . 8929E-04 
0 . 8997E-04 
0 . 9065E-04 
0.9132E-04 
0.9199E-04 
0 . 926SE-04 
0 . 9332E-04 
0 . 9398E-04 
0 . 9483E-04 
0 . 9529E-04 
0 . 9594E-04 
0 . 9858E-04 
0 . 9723E-04 
0 . 9788E-04 
0 . 9850B-04 
0.9914E-04 


0 . 4935B-04 
0.5038B~04 
0 . 5141E-04 
0 . 5244E-04 
0 . 5347E-04 
0 . 5450B-04 
0 . 5553E-04 
0 . 5868E-04 
0 . 5780E-04 
0 . 5883E-04 
0 . 5988B-04 
0 . 8070B-04 
0 . 8173B-04 
0 . 8278E-04 
0 . 8380E-04 
0 . 8484E-'04 
0 . 8S87E-04 
0 . 8891E-04 
0.8794E-04 
0 . 8898E-04 
0.7001E-04 
0.7105E-04 
b . 7209E-04 
0 . 7313E-04 
0.7418E-04 
0.7520E-04 
0.7824E-04 
0 . 7728E-04 
0.7831E-04 
0.793SB-04 
0.8039E-04 
0.8143E-04 
0.8247E-04 
0.8351E-04 
0 . 8455E-04 
0 . 8559E-04 
0 . 8883E-04 
0 . 8787E-04 
0.8871E-04 
0 . 8975E-04 
0.9079E-04 
0.9183B-04 
0 . 9287E-04 
0.9391E-04 
0 . 9495E-04 
0.9599E-04 
0.9703E-04 
0 . 9808E-04 
0.9912E-04 
0 . 1002E-03 
0.1012B-03 
0 . 1022E-03 
0 . 1033E-03 
0 . 1043E-03 
0 . 1054E-03 
0 . 1084E-03 
0 . 1075E-03 
0 . 1085E-03 
0 . 1095E-03 
0.1108E-03 


0.6650B-t>01 

0.6600E-»-01 

0.6650B-»^01 

0.6700B-»-01 

0.6750B-t>01 

0.6800B-«>01 

0.6850B-«-01 

0.6000B-»-01 

o.egsoB+oi 

0.7000B+01 

0.7050B+01 

0.7100B+01 

0.71S0B-t-01 

0.7200B+01 

0.72S0B+01 

0.7300B+01 

0.7350B+01. 

0.7400B+01 

0.7450B+01 

0.7500B+01 

0.7550B+01 

0.7600E+01 

0.7650B+01 

0.7700E+01 

0.7750B+01 

0.7800E+01 

0.7850E+01 

0.7900E+01 

0.7950E+01 

0.8000B+01 


0.0076B-04 
0.1004E-03 
0 . 1010E>03 
0 . 1016E-03 
0 . 1022E-03 
0 . 10298-03 
0 . 1035E-03 
0.1041E-03 
0 . 1047E-03 
0 . 1053E-03 
0 . 1059E-03 
0.1065E-03 
0 . 1071E-03 
0 . 1077B-03 
0 . 1083E-03 
0 . 1088E-03 
0 . 1094E-03 
O.llOOB-03 
0.1106E-03 
O.llllE-03 
0.1117E-03 
0.1123E-03 
0.1129E-03 
0.1134E-03 
0.1140E-03 
0.1145B-03 
0,1151E-03 
0.1156E-03 
0.1162E-03 
0.1167B-03 


0.1116E-03 
0.1127E-03 
0.1137E-03 
0.1148E-03 
0.11S8E-03 
0.1168E-03 
0.1179B-03 
0.1189E-03 
0 . 1200E-03 
0 . 1210E-03 
0 . 1221E-03 
0 . 1231E-03 
0 . 1242E-03 
0 . 1252E-03 
0 . 1262E-03 
0 . 1273E-03 
0 . 1283E-03 
0 . 1294E-03 
0 . 1304E-03 
0 . 131SE-03 
0 . 1325E-03 
0 . 1336E-03 
0 . 1346E-03 
0 . 13S7E-03 
0 . 1367E-03 
0 . 1378E-03 
0 . 1388E-03 
0 . 13g8E-03 
0 . 1409E-03 
0.1419E-03 


c  This  main  program  executes  a  Curve  Fit  using: 
c 

c  STEPIT  for  the  curve  fitting 

c  STSBT  for  initialization  of  STEPIT  input  parameters 

c  MODEL  for  minimizing  the  current  equation  derived 

c  from  a  "Charge-Sheet  Model" 

c  CSO  for  main  Charge-sheet  program 

c  FUNCTION  for  the  current  equation 

c 
c 

c  The  user  must  supply: 

c 

c  1.  IV  data  file  in  required  format 

c 

c  2.  Parameter  file  in  required  format 

c 

c  3.  Initial  Starting  value  file  for  mobility 

c  and  flat  band  voltage  in  each  data  set. 

c 
c 

C 

c  The  user  can  request  a  copy  of  STEPIT  from: 

c 

c  Quant\im  Chemistry  Program  Exchange 

c  Indiana  University 

c  Bloomington,  Indiana  47401 

c 

c  Authored  by:  J.P.  Chandler 

c  Department  of  Computing  and  Information  Sciences 

c  Oklahoma  State  University 

c  Stillwater,  Oklaihoma  74074 

c 
c 

C  *  I** ^  ’t*  4> 4* ^ 4' 4*  ■t'*’** Ml  ^  *  :4i  :4i  :4i  ^ 

C 

c 

c  The  program  CSFIT  assumes  that  the  fitting  parameters 

c  are  flat  band  voltage  and  mobility  with  the  flat  bamd 

c  voltage=X(2i-l)  and  mobility=X(2i)  where  iM:l,2, . 

c  .  .  .nvds=the  total  niuaber  of  data  sets  and  2Minvds=total 

c  number  of  fitting  parameters.  Thus,  if  2  I-V  data  sets 

c  a.re  used  with  X(l)=ul,  X(2)=fbvl,  X(3)=u2,  and  X(4)=fbv2 

c  the  output  file  results  would  be  as  follows: 

c 
c 

clSUBROUTINB  STEPIT.  COPYRIGHT  (C)  1965  J.  P.  CHANDLER 
c 

c  INITIAL  VALUES .... 


c 

c 


c 

MASK 

= 

0 

0 

0 

0 

c 

c 

X 

0 . lOOOOE+00 

0 . 50000B+03 

O.lOOOOE+00 

0 . 50000E+03 

c 

c 

XMAX 

= 

0 . lOOOOE+36 

0 . lOOOOE+36 

0 . lOOOOE+36 

0 . lOOOOE+36 

c 

c 

XMIN 

s 

-0 . lOOOOE+36 

-0 . lOOOOE+36 

-O.lOOOOE+36 

-0 . lOOOOE+36 

c 

24  c 

DELTX 

— 

0 . lOOOOE-02 

0 . 50000E+01 

O.lOOOOE-02 

0.50000E+01 

DELMN  =  O.lOOOOE-03  0.10000E--03  O.lOOOOE-03  O.lOOOOE-03 


4  VARIABLES,  4  ACTIVE.  MATRX  =  110 

NFLAT  =  1  RELAC  =  O.lOOOE-16 


NFMAX  =  1000000 


FOBJ  =  0.5310057276B-08 


BEGIN  MINIMIZATION. 


TERMINATED  WHEN  THE  STEP  SIZES  BECAME  AS  SMALL  AS  THE  DELMN (J) . 


277  FUNCTION  COMPUTATIONS 


FINAL  VALUE  OF  FOBJ  =  0 . 136086074S97773B-10 


FINAL  VALUES  OF  X( J) .... 


0.897636003046762E+00  0 . S38229046900407E+03  0.903491922574461E+00 

0 . 537289871078474E+03 

ISUBROUTINE  STERR. 

COMPUTE  AN  APPROXIMATE  ERROR  MATRIX  USING  FINITE  DIFFERENCES. 

INCREMENTS  IN  X( J)  TO  BE  USED .... 

0 . 32000E-t-02  0.20000E+05  0.32000E+02  0 .  lOOOOE+OS 

SQUARE  ROOT  OF  NEGATIVE  NUMBER  HAS  OCCURRED: 

EITHER  THE  INITIAL  STARTING  VALUES  ARE  BAD  AND  NEED  TO  BE  RESET  OR 
THE  ERROR  MATRIX  CAN  NOT  BE  COMPUTED  WITH  THE  CSO  ALGORITHM, 

THUS  IGNORE  MESSAGE  AND  ASSUME  FITTED  VALUES  ARE  CORRECT. 


Note:  The  error  matrix  could  not  be  computed  and  thus  it  is 
disregarded.  When  an  approximate  error  matrix  is 
attainable,  a  matrix  of  the  second  partial  derivatives, 
a  standard  errors  matrix,  and  a  lower  triangle  of  the 
correlation  matrix  is  given. 
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I  '  1^  I 


*llR' 


p 


c 

c 


The  final  optimized  values  are  given  in  corresponding  order, 
vfb(l)  ,mob(l)  ,vfb(2)  ,mob(2)  . . vfb(i)  ,fflob(i)  . 


i.e. 


% 

I 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


J 

> 


200 

c 


implicit  realms  (a-h,o-z) 
external  model 
dimension  a (1000, 3) 

common  /cstep/  x(20) ,xmax(20) ,xmin(20) ,deltx(20) ,delmn(20) , 
*err(20,21)  ,fobj  ,nv,ntrac,matrx,mask(20)  ,  nf  max,  n^at,  j  vary , 
♦nxtr a , kf lag , nor ep , ker f 1 , kw 
common  /cmod/  xm(ioOO) ,y(1000) ,npts 
common  /input/  npset (20) , vds (20) ,nvds 
real>«i4  vds, beta 
common  /mshape/  djxin, ljun, Ich 
real>ti4  djun,  1  jun, Ich 
common  /dopant/  nsd,nch 
real *4  nsd,nch 

common  /bias/  vd,vs,vg,vsub,vfb 
real *4  vd,vs,vg,vsub,vfb 
common  /gate/  tox,q0,kox,k8i 
real>t>4  tox,q0,kox,ksi 

common  /sheet/  mustar , z,p8ib, Ib,cox,p8is0,p8isl ,n0, iO, id 

real>t>4  mustar ,  lb , nO ,  iO ,  id , psib , cox, psisO , psisl ,  z 

COMMON  /NORM/  VT,VNORM,Q,NI,EPSO,BPS,LNORM,NNORM, tnorm,muO 

REAL>o4  VT ,  VNORM ,  Q ,  NX ,  EPSO ,  EPS ,  LNORM ,  NNORM ,  tnorm ,  muO 

character'«i40  f ilinl ,f ilin2,f iiin3,filout 

write  (6,*)  ’Enter  IV  Data  file* 

read  (5,5)  filinl 

format  (a) 

write  (6,*)  ’Enter  Parameter  file’ 
read  (5,5)  filin2 

write  (6,>t>)  ’Enter  Initial  Starting  Value  file’ 

read  (5,5)  filinS 

write  (O,*)  ’Enter  Output  file’ 

read  (5,5)  filout 

open  (unit=10 , f ile^f ilinl , status^ ’ old ’ ) 
open  (unit=ll , f ile-f ilin2 , status^ ’old ’ ) 
open  (unit=12, f ilesf ilin3, status=’old ’ ) 
open  (imitslO ,  f ile=f ilout ,  status^  ’  new  ’ ) 
open  (uni t=l7,file=’ plot.dat’ , status= ’ new ’ ) 

Read  from  the  Data  file.  The  following 
format  for  the  data  file  is  required: 


npts 

a(j  »1) 

where : 


a(j,2) 


a(j,3) 


nptsstotal  number  of  iv  points 
a(j > l)=source-to-drain  voltage 
a(j,2)=gate  voltage 
a(j)3)=drain  current 


read(10,>«>)  npts 
do  200  3 si, npts 

read  (10, ♦)  a(j  ,  1)  , a(j  ,2)  ,a(j  ,3) 
y(j)=a(j  ,3) 

3tm(j)=a(3  ,2) 
continue 


/.  V. 


V-V .JW  ^  >  ^ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

201 

202 


c 

c 

c 

c 

c 

c 

c 

c 


Read  from  the  Parameter  file.  The  following 
format  is  required: 

nwds 

▼ds(l) 

▼d8(i) 

npset(l) 


npset (i) 

nch 

tox 

1  jun 

Ich 

z 

where :  nvds=the  number  of  data  sets . 

vds (i)=source-to-drain  voltages  for 
each  set. 

npsets (i)=nufflber  of  data  points  in 
each  set. 

nch=avera,ge  channel  doping  (cm-S) 
tox=oxide  thickness  (cm) 

Ij unilateral  source  and  drain 
diffusion  (cm) 

Ichidrawn  length  of  polysilicon  (cm) 
*  Ziwidth  of  channel  (cm) 

read  (11,*)  nvds 

do  201  j lil , nvds 

read  (11,*)  vd8(jl) 

do  202  j2il,nvds 

read  (11,*)  npset(j2) 

read  (11,*)  nch, tox, 1 jun, Ich, z 

q=1.6E-19 

ni=l  .5E-«-10 

nchinch/ni 

koxiS . 9 

ksisll .8 

eps0=8 . 854E-14 

muOilOOO 

q0=0. 

vt= . 0259 

betail .0/vt 

Inormisqrt (ksi*eps0/q/ (beta*ni)) 

toxi tox/ Inorm 

1 j un= 1 j un/ Inorm 

Ichilch/lnorm 

Ziz/lnorm 


Call  STSET  to  initialize  starting  values  for  STEPIT. 
See  STEPIT  documentation  for  details. 

Note:  STEPIT  allows  the  user  to  control  the  least 

squares  criterion  by  setting  certain  values: 

x(i)iinitial  values  of  the  parameters,  on  which 
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the  function  to  be  minimized  depends . 
xmajc(i)3:the  upper  limit  on  the  value  of  x(i)  . 
xmin(i)sthe  lower  limit  on  the  value  of  x(i) . 
deltx(i)sthe  initial  step  by  which  x(i)  is  varied. 
delmn(i)s:the  lower  limit  (convergence  tolerance) 
on  the  step  size  for  x(i) . 
mask(i)=set  nonzero  if  x(i)  is  not  to  be  varied, 
but  rather  is  to  be  held  at  its  initial 
value . 

ntrac^determines  the  amount  of  printout  from  STEPIT 

ntracsl  to  obtain  trace  during  the 
minimization  process 

ntracsO  for  initial  and  file  output  only 

ntracs-1  for  no  output  except  error  messages 

STSET  will  set  all  the  above  to  default  values  except  ntrac 

After  SESET  is  called  the  inital  starting  values  are  read 
in  via  the  following  format. 

initial  vfb(l) 
initial  mob(l) 
initial  vfb(2) 
initial  mob (2) 


initial  vfb(nvds) 
initial  mob(nvds) 

call  stset 

nv=nvds*2 

ntracsO 

do  203  j3sl,nv,2 
read(12,«)  x(j3)  ,x(j3+l) 
x(j3+l)  =  -x(j3+l) 
do  204  j  4=1 , nv 
delmn ( j  4) =1 . OE-04 
matrx=110 

call  STEPIT  for  least  squares  fit. 
call  stepit (model) 


print  results 


k=l 

do  2000  j=l,nvds 

Vffit  =  -x(2*j-l) 

write (16, 1002)  Vff it,x(2*j) ,vds(j) , j 

format  (’  Vfb  »  ’,lpel2.6,’  Mu  -  ’,0pf7.2,’  Vd  =  ’,f7.2, 
’  Set  ’,i2,/,’  Vg  Id  Ideal’) 

write (17,*)  npset(j) 
imax=npset  (  j  ) 
vd=vds(j)/vt 


▼£b»-x(2*3“l)  /▼■t 

muatarsx (2* j ) /muO 
do  1000  i=l,imax 
vg=oaii(k) /vt 
call  csO 

write (16, 1001)  xm(k) ,y(k) ,id 
1001  format  (3el2.4) 

write (17, 1001)  ▼ds(j) ,xm(k) ,id 
ksk+1 

1000  continue 

2000  continue 

207  continue 
close (10) 
close (11) 
close (12) 
close (16) 
end 
c 
c 

c  MODEL  subroutine 

c 

c 

subroutine  model 
implicit  realms  (a-h,o-z) 

common  /cstep/  x(20) ,xmax(20) ,xmin(20) ,deltx(20) ,delffln(20) , 
*err  (20,21)  ,fobj  ,nv,ntrac,matrx,fflask(20)  ,nfmax,nf  lat,  j  vairy , 
•^nxtr  a ,  kf  lag ,  norep ,  ker  f  1 ,  kw 
common  /cmod/  xm (1000) ,y (1000) ,npts 
common  /input/  npset (20) ,vds (20) ,nvds 
real*4  yds, beta 
common  /mshape/  dj\m,ljun,lch. 
real *4  djun,  1  jxin, Ich 
COMMON  /GATE/  TOX,QO,KOX,KSI 
REAL*4  K0X,KSI,tox,q0 
common  /dopant/  nsd,nch 
real >•>4  nsd,nch 

COMMON  /NORM/  VT,VNORM,q,NI,BPSO.EPS,LNORM,NNORM, tnorm,mu0 
RRAL*4  VT,VN0RM,q,NI,EPS0,EPS,LNbRM,NN0RM,tnorm,mu0 
common  /bias/  vd,vs,vg,vsub,yfb 
real *4  vd,ys,vg,vsub,vfb 

coimnon  /sheet/  mustar,z,psib,lb,cox,psis0,psisl,n0, iO, id 

real4i4  mustar ,  lb,n0,  iO,  id,psib, cox,psis0, psisl ,  z 

real>ti4  if\in,ni 

fobj=0. 

k=l 

do  2000  jsl,nyds 
imaxsnpse t ( j ) 
vd=vds(j)/yt 
vfb=-x(2#j-l)/vt 
mustarsx(2* j ) /muO 
do  1000  i=l , imax 
▼g*xm(k)/vt 
call  csO 

fob  j  =f  ob  j  +  (id/y  (k)  -1 . 0)  *■•2 
k»k-t>l 
continue 
continue 
return 
end 


1000 

2000 


CSO  Subroutine 


subroutine  csO 

common  /msbape/  djun, lj\in,lcb 
reml*4  djtin,  1  jun, Icb 
COUMON  /GATE/  TOX,QO,KOX,KSI 
REAL4i4  KOX,KSI,tox,qO 
common  /dopant/  nsd , neb 
real4>4  nsd, neb 

COMMON  /NORM/  VT , VNORM , Q , NI , EPSO , EPS , LNORM , NNORM , tnonn , muO 
REAL*4  VT, VNORM, Q, NI, EPSO, EPS, LNORM, NNORM, tnorm.muO 
common  /bias/  vd,TS, vg, vsub,Tfb 
real *4  vd,vs, vg,vaub,vfb 

common  /sheet/  mustar,s,psib,lb,cox,psisO,psisl,nO,iO,id 

real *4  mustar , lb,nO, iO, id,psib, cox,psisO,psisl , z 

real*4  ifun,nl 

beta  *  1 . 0/vt 

psib  -  ▼t’t'alog  (neb) 

lb  =  sqrt(ksi4>eps0/q/(beta*ncb>*>ni)) 
cox  3  kox4ieps0/ (T0X*LiN0RM) 
vgl  =  vg-vfb+q*beta*QO/cox 

bpsO  =  2 .0*beta4ipaib-2 .0>tialog(beta>osqrt(2 .0)  *q>t>ncb*ni>*ilb/cox) 
bpsi  3  bpsO 
bpsis  3  l.OelO 

if  (vgl . It .bpsO)  bpsis  3  bpsO 
if  (vgl . le .bpsO)  goto  11 
do  10  i»l,500 

bpsi  3  bpa0+alog(abs((vgl-bp8i)*32/(l.+(bpsi-1.0)/exp(bp8i 
1  -2.0«ip8ibi«beta)))) 

if  (abs(bp8i-bpsi8) . It . 1 .e-4)  goto  11 
bpsis  3  bpsi 
continue 
psi803bpai/beta 

nOsncb^ni^ilb^sqrt  (2  . )  ■*  (  (bpei-1 .0+  (1 .  /neb)  *>*2* 

1  exp (bpsi) ) *30 . 5- (bpsi-l . ) **0 . 5) 

psisl  3  p8is0-t-vd*vt 
vO  3  (q*ncb*ni*lb/cox) 3*2 
IP  (VGl.LT.  1.0)  GO  TO  5000 

psiss3vgl3vt-fv03beta-8qrt  (  (v03beta)  3*2+2 . 3v0*  (vgl-1 . 0)  ) 
if  (psisl .gt.psiss)  psislspsiss 
nl  3  0  J 

if  (ps...sl .  eq. psiss)  goto  13 
b  *  1.0 
dpold  3  1.0 
do  12  i3l,lOO 

nl  3  cox3(vgl3vt-psisl) /q-ni3ncb3lb38qrt (2 . ) * (beta3psi8l-l . ) **0 , 
if  (nl.lt.l.Oe-25)  nl*1.0e-25 
psisi3psia0+vd3vt+alog (nl/nO) /beta 
delp8i3p8isi -psisl 

if  (ab8(delp8i) .gt.vt)  delp8i3vt*ab8 (delpsi) /delpsi 

if  (delpsi . eq. 0 .0)  goto  12 

if  (dpold/delpsi . It .0.0)  h3b/2.0 

psisl3psisl+b3delp8i 

if  (abs (delpsi) . It. 1 .Oe-5)  goto  13 

dpold3delp8i 

continue 

id  3  ifun (psisl , Icb) 


xidot  s  id>tibeta4i>t>2/(mu8tax<t>mu0«z/(lcli-2 j\m) -fcox) 
xlch  s  lcli-2.0>«iljTin 

c  WRITE  (6,90)  BBTA,PSIB,LB,C0X,BPS0,BPSI,PSIS0,N0 
c  1  , PSISL,NL, ID, iO, xidot, vg,Yd, xlch 

90  FORMAT  (’  BETA  =  MPE15.6/ 


1 

) 

PSIB 

S 

’,1PB15.6/ 

1 

9 

LB 

s 

’,1PB15.6/ 

1 

9 

COX 

s 

’,1PE15.6/ 

1 

9 

BPSO 

s 

’ , 1PB15 . 6/ 

1 

9 

BPSI 

s 

’  ,1PB15.6/ 

1 

9 

PS(0) 

= 

’ , 1PB15 . 6/ 

1 

9 

NO 

s 

’ ,1PB15.6/ 

1 

9 

PS(L) 

= 

’ , 1PB15 . 6/ 

1 

9 

NL 

= 

’  ,1PE15.6/ 

1 

9 

ID 

= 

’,1PB1S.6/ 

1 

9 

10 

’ , 1PE15 . 6/ 

1 

9 

I.  .  . 

’ ,1PE15.6/ 

1 

9 

VG 

’ ,1PE15.6/ 

1 

9 

VSD 

= 

’ , 1PE15 . 6/ 

1 

9 

LCH 

’ , 1PB15 . 6/) 

return 


5000  id  =  l.Oe-20 
cSOOO  write (16, 500) 

500  FORMAT (//, ’SQUARE  ROOT  OF  NEGATTVE  NUMBER  HAS’, 
1’  OCCURRED:  ’//’EITHER  THE  INITIAL  STARTING’, 

1 ’  VALUES  ARE  BAD  AND  NEED  TO  BE  RESET  OR ’ , 

1//, ’THE  ERROR  MATRIX  CAN  NOT  BE  COMPUTED  WITH’ , 
1’  THE  CSO  ALGORITHM, ’//’THUS  IGNORE  MESSAGE’, 

1’  AND  ASSUME  FITTED  VALUES  ARE  CORRECT. ’) 
c  STOP 

return 
end 


c  Function  Subroutine  IFUN 


real  function  ifun(ps,y) 
common  /mshape /  d j  un , 1 j  un , Icb 
real 4^4  djun,  ljun,lcb 
COMMON  /GATE/  TOX,qO,KOX,KSI 
REALti4  KOX,KSI,tox,qO 
common  /dopant/  nsd,ncb 
real *4  nsd,ncb 

COMMON  /NORM/  VT,VNORM,Q,NI,EPSO,EPS,LNORM,NNORM, tnorm,muO 
REAL>*«4  VT ,  VNORM ,  Q ,  NI ,  BPSO ,  EPS ,  LNORM ,  NNORM ,  tnorm ,  muO 
common  /biaa/  vd, vs ,Yg, vsub , Yfb 
real*4  vd,vs, vg,vsub, vfb 

common  /sheet/  mustar , z , psib, lb , cox, psisO, psisl ,nO, iO, id 
real*4  mustar , lb , nO , iO , id , z , psib , cox , psisO , psisl 
real ■<•4  iOO,fff 
beta  =  1.0/vt 

vgl  =  vg-vfb+q>t«beta*QO/cox 

((y-2.0>*«lj\in)  .ne  .0.0)  iOO  =  1 .0/beta*mu8tar>*'mu0x«z/(y-2 .0*1  jun) 
fff  *  cox* (1 .0+vgl) * (ps-psisO) -beta/2 . *cox* (p8**2-p8is0**2) 
f £f =f£f-q*ni*NCH*lb*8qrt (8 .) /3. * ( (beta*ps-l .)**!. 5-(beta*PSIS0 
1  -l.)**1.5) 

£££*££ £+q*nch*ni*lb*8qrt (2 . ) * ( (beta*ps-l . ) ** . 5- (beta*psisO-l . ) 

1  **.5) 

i£un  ■  iOO*£££ 
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